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We develop a semiclassical theory of the nondegenerate parametric amplification in 
a single miniband of superlattice. We present the formulas describing absorption and 
gain of signal and idler fields in superlattice and analyze the limiting cases of strong 
and weak dissipation. We show how the well-known Manley-Rowe relations arise in the 
tight-binding lattice in the weak dissipation limit. Our results can be applied to an 
amplification of THz signals in semiconductor superlattices and a control of nonlinear 
transport of cold atoms in optical lattices. 
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1. Introduction 

Terahertz radiation (0.3 — 10 THz) has enormous number of promising applications 
in different areas of science and technology, not least because certain important 
molecular resonances occur at THz frequencies. One of the main challenges is to 
construct a coherent miniature source/ amplifier of THz radiation that can operate 
at room temperature. 

The classical works of Esaki and TsrP and Ktitorov et al® stimulated enormous 
activity devoted to the non-linear interaction of a high-frequency electric field with 
miniband electrons in dc-biased superlattice. Bloch oscillations of miniband elec- 
trons in a such kind of system can, in principle, result in an amplification of THz 
radiation (for a review, see Ref. [3]). The main problem in the realization of these 
THz amplifiers is a formation of high-field electric dom ains inside a superlattice 
under conditions of negative differential conduct i vity^^ (NDC). Electric domains 
are destructive for gain at THz frequencies.^ Several novel ideas have been recently 
introduced in order to solve the problem (see Ref. [7] and papers cited therein). 



November 
proc'L3'corr'eprint 



10, 



2008 11:52 WSPC/INSTRUCTION 



FILE 
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Another suggested type of THz amplifiers with a superlattice as an active ele- 
ment is based on a coherent interaction of alternating elec tric fields with two com- 
mensurate frequencies (Refs. IMTB]). In recent Letters™ we clarified the paramet- 
ric nature of this effect and found that it can exist without switching the operation 
point to the NDC part of the time-average voltage-current characteristic of the 
superlattice. This theoretical result allows to expect that the undesirable electric 
instability can be effectively suppressed in the case of parametric gain. 

The parametric amplifiers are characterized by a time-variation of a reactive 
parameter which pumps energy into an amplified signal, i f th e phases and the fre- 
quencies of the signal and parameter are suitably related!^ In a single miniband 
superlattice the para met ric resonance occurs due to oscillations of the effective mass 
of miniband electrons^, which is caused by the pump field 



In such degenerate parametric amplifier there is only one resonator mode whose 
frequency io± is a harmonic or half harmonic of ui. The degenerate amplification in 
superlatti ces has be en analyzed in several physically interesting limits, including the 
quasistatiJUE^Hl 

(ujt <C 1, WiT <C 1) and "semiquasistatic" (ujt < 1, 

u)\T > 1, r is the intraband relaxation time). 

The degenerate parametric amplification is a phase-sensitive process. On the 
other hand, nondegenerate parametric amplifiers can provide a phase-insensitive 
amplification of a given signal!^ In this case there are two resonator modes with 
incommensurate frequencies ui\ and uii satisfying u>\ ± u>2 = nui (n is an integer). 
There still exist a few p ublica tions devoted to the nondegenerate parametric ampli- 
fication in superlattices 1181333 

In this paper, we re-examine the problem of nondegenerate parametric amplifi- 
cation in single miniband of superlattices and analyze the limits of weak (r — > oo) 
and strong (r — > 0) dissipation. In the limit of weak dissipation, we show that our 
formulas for absorption and gain satisfy the well-known Manley-Rowe relationsPSl 
We concentrate mainly on the small-signal amplification and demonstrate that the 
net absorptions of weak probe fields in superlattice always consist of two distinct 
terms describing the incoherent and parametric mechanisms. We find that for both 
frequency relations u>i ± u>2 = nu> and odd n the parametric term is zero in unbi- 
ased (i?dc = 0) superlattice. Our results provide a mathematical background also 
for analysis of the nondegenerate parametric amplification in various other physical 
systems demonstrating nonlinear transport in a single energy band. 
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Nondegenerate parametric amplification in superlattices 3 

2. Main formulas describing the nondegenerate parametric 
amplification in a superlattice 

The stationary current (t ^> r) in the tight-binding miniband can be calculated for 
an arbitrary time-dependence of the electric field E(t) aJ^ 



j{t) = 2i peak Im J F(i), 



F(t) = i f dt 1 exp^-(t-t 1 )/T + l I'dh^^ | 



(2) 



where j pea k is the peak current density in the Esaki-Tsu current-field characteristic 

NeAdh(A/2k B T) 
Jpcak Ah I (A/2k B TY [ ' 

A is the miniband width, d is the period of the superlattice potential, N is the 
density of carriers in the first superlattice miniband, I n {x) are the modified Bessel 
functions and T is the lattice temperature. 

In our case the total electric field acting on miniband electrons is 



E(t) = E pnmp (t) + Ei cos (wit + <f>i) + E 2 cos (cj 2 t + <f> 2 ). 
By straightforward calculation we get 

F= J2 HP)Ji+h{P)Jk{Pi)Jk+c{Px)J m {fo)J m+ 02) 

l,k,m,h,c,j 



(4) 



exp 



ihtot + ic{u)\t + 4>\) + ij(u)2t + 4> 2 ) 



1 — iuJsT — ilu>T — ikuoiT ~ imui 2 T 
where Ji(x) are the Bessel functions, /3 = eEod/hiu, = eEid/hwi (i — 1,2) and 
u>b = sEdcd/h is the Bloch frequency. All summations are from — oo to oo. 
Thus 

j(t) = MP)Ji+h(P)Jk(Pi)Jk+c(Pi)J m (l32)J m +M 

l : k 7 m,h,c 7 j 

x jet(^BT + loJT + koJ\T + mu>2T) COS (htot + c(u}%t + 4>\) + j(u>2t + 4>2)) 

+K(u>bt + l^T + kuj\T + muj 2 T) sin [hut + c(uj\t + 4>i) + j{^it + (f> 2 )) 



(5) 



where j et is the Esaki-Tsu current-field characteristic 

. , v 2ui b t 

Jet{VBT) = Jpeak— — j yy 

1 + {lo b tY 

and K is connected to the j e t by the Kramcrs-Kronig relations and is given by 

2 



(6) 



K(u B t) = jpcak 



+ (uBTf 



(7) 
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The formula §5§ describes the time-dependent current response in the case of trichro- 
matic electric held, and both the dc current jd c — (j(t))t and the absorptions at 
different frequencies A u = (j(t) coscut) t and A Wi = (j(t) cos((Jit+<fii))t {i = 1, 2) can 
be calculated from this formula. Here (. . .)( denotes the time-average. The gener- 
ated/absorbed power densities at different frequencies are then given by = AujEq 
and V UJi — AujtEi (i — 1,2). Absorption and emission correspond to V > and 
V < 0, respectively. 



2.1. Absorptions for arbitrary amplitudes of the ac fields 
satisfying u> 1 + uj 2 — nui 

By assuming that the mutually incommensurate frequencies satisfy the condition of 
parametric resonance lj± + lu 2 = nuj, we get from formula ([5]) the following formulas 
for the absorptions at different frequencies 



COS (j(>l + 02 )) 



+K(lubt + Ilut + kujiT + miti-iT) 



sin (j((j>i + <f> 2 )) 



(8) 



A^ = JiWVi-MJM 



j et {0JBT + lWT + klO\T + mU>2T) 



cos (j(0i +(j> 2 )) 



Jm{(32)Jm+j(P2) 



-K(lubt + Ilut + koj\T + mw2i")- 



si.n + <j> 2 )) 



(9) 



and 



l,j : k : m 

COS (j(0i + 2 )) 



Jm+j+1 {P2) + Jm+j-lifo) 



jet(u>BT + lVT + kuJlT + mUJ 2 T) 



-K{(jJbT + lljJT + k<jJ\T + muj 2 T) 



sm 



(10) 
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Nondegenerate parametric amplification in superlattices 5 

2.2. Absorptions for arbitrary amplitudes of the ac fields 
satisfying Wi — lj 2 — noo 

Similarly one gets in the case of lo\ — lo 2 = nui the following formulas for the 
absorptions at different frequencies 



l,j,k,m 



Jl-jn+l(P) + Jl-jn-l(P) 



Jk (Pi ) Jk+j (Pi ) Jm (P2 ) Jm-j (Pi ) 



j et (u>BT + Ilot + kuiT + moo 2 T 



+K(ojbt + Ilot + kuo\T + mL0 2 T) 



cos - <l>2)) 



sin - <j) 2 )) 



(11) 



•At, 



and 



]T Jl(P)Jl-jn(P)Jk(Pl) 
l,j,k,m 



jet(iOBT + IlOT + kiJiT + mW 2 T) 



Jk+j+l(Pl) + Jk+j-l(Pl) 
COS - 2 )) 



J m (P2)J m -j(p2) 



+K(lu b t + Ilot + kuoiT + rriL0 2 T) 



sin (j'Oi - fa)) 



(12) 



Y, MP) Jl-jn (P) Jk (Pi ) Jk+j (Pi )Jm(p2) 

l,j,k : m 

cos (jOi - (j} 2 )) 



Jm-j+l(P2) + J m -j-l(p2) 



jet(iOBT + Ilot + kiO\T + mu> 2 T)- 



+K(lobt + Ilot + ku\T + itiuo 2 t) 



sin(j(>i -<f) 2 )) 



(13) 



2.3. Manley-Rowe relations in the limit of weak dissipation 

In conditions of parametric resonance lo\ + lo 2 = noo in unbiased superlattice, we 
get in the limit of weak dissipation r — > 00 



k,j 



J-kn-jn+l(P) + J-kn-jn-l(P) 
XjpoakSin {j((j)l +(fa)), 



Jk(Pl)Jk +j (Pl)Jk(p2)Jk +j (p2) 

(14) 



A 0Jl = ^ J-kn(P)J-kn-jn(P)Jk(Pl 
k,j 

XjpoakSin (j((j>l + (j} 2 )) 



Jk+j+l(pl) + Jk+j-l(pl] 



Jk(p2)Jk+j(p2) 

(15) 
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and 

kn—jn 

xjpoaksin (j(0i + h))- (16) 

By using the Bessel function relation J„_i(a;) + J n+1 {x) = 2nJ n (x)/x, we get the 
following relations 

Vu/nw = -V Ul /wi, 
Vujjnuj = —Voj 2 /w2, 

V Ul /ui=V Ua /w2, (17) 

where the last equation actually follows from the first two relations. These sim- 
ple connections between the generated/absorbed power at different frequencies are 
known as the Manley-Rowe relations and are typical for parametric systems consist- 
ing of nonlinear inductors and capacitors^]. Indeed in the limit of weak dissipation, 
electrons in superlattices are gathered together in the quasimomentum space form- 
ing an electron bunch which oscillates coherently in quasimomentum space along a 

ipi I 

certain ballistic trajectory^ 1 . Here the effective mass of the electrons in the bunch 
plays the role of the differential inductance leading to phase-dependent small-signal 
absorption and gain for the probe fields if the condition of parametric resonance is 
satisfied. 

Similarly in conditions of parametric resonance lo\ — u>2 = nuj, the Manley-Rowe 
relations in unbiased superlattice are given by equations 

P^/nuj = -Vu,Ju\ 

Vujui = -V u Ju 2 . (18) 

3. Absorption of weak signal and idler fields in the presence of 
strong pump 

If the pump is strong but the probe fields (signal and idler) are weak, modifications 
of the voltage-current characteristic of superlattice are determined by the action 
of the pump field alone. In this case, the problem of electric stabilit y has been 
considered earlier in the theory of degenerate parametric amplificationff2ll2l 

Now we want to find the expressions for absorptions of weak probe fields <C 1 
(i = 1,2). By using the Bessel function expansion [Jq(x) rs 1, J±i(x) « ±x/2], we 
find from the formulas ([9]), (fl7J|) . (fT2|) and (fl3|) that absorptions at frequencies uj\ 
and Co>2 can be written as 

a yiincoh _i_ jcoh 

A^=AZT h + AZf. (19) 
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3.1. Incoherent components of absorption 

The phase-independent incoherent absorptions are always given by 



j et {iOBT + IUT + LUit) - j e t{^BT + lUT - UJ\t) 
jet(w B T + l0JT + W 2 t) - jeti^BT + l^T - U 2 t) 



(20) 



Without the pump these equations describe a free-electron absorption in unbiased 
superlattices-^. On the other hand in the presence of the pump field, they can in 



tne pump t 
ticesffl On 

principle also describe the Bloch gairPl 

The phase-dependent coherent absorptions result from the parametric resonance 
in superlattice minibancP^ and depend on the relation of the frequencies. 

3.2. Coherent absorptions for u>± + u>2 = w 

For uj\ + u>2 = nu the coherent part of the absorptions can be written as 
-Cf = cos(0i + fa) + *C Mfa + fa) 



Alt = B™cos(fa 



B s £sin(fa + fa), 



where 



K7COS 
D U1 1 



Y,HP)Jl-n{P) 



jet{^BT + iLOT + UJ 2 t) - jeti^BT + liOT 
jeti^BT + liOT) - j et {uJBT + liOT - LU 2 t) 



K(u>bt + Iu>t) — K(lobt + luT + lu 2 t) 
K(ujbt + Ilot) — K{lobt + £wt — lu 2 t) 



i^Ji(/3)J ; _ n (/3) 



jeti^BT + lUT + Wit) - j et (w B T + IuJt) 
jet{u B T + Iu>t) - jeti^BT + liOT - W].t) 



and 



K{ujbt + Ilot) — K(lobt + iwr + wit) 
K(ujbt + Zwt) — K{lobt + iwr — Wit) 



(21) 



(22) 



(23) 



(24) 



(25) 
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3.3. Coherent absorptions for Wi — u>2 = nw 

On the other hand for uj\ — u>2 = nuj the coherent part of the absorptions can be 
written as 



/JCOh 



Alt = »£ s cos(^2 -<h) + M<h -(h), 



(26) 



where 



+xE J 'C 8 ) J i-nC8) 



j et (uJBT + ILOT) - j et (uJBT + IUJT - UJ 2 t) 
JetiuiBT + ILOT + LU 2 t) - jet&BT + IlOt) 



(27) 



B sin 



^J,(/3)J i+n (/3) 



K(ujbt + Iujt) — K{u>bt + Iujt — UJ 2 t) 
K{u>bt + Iujt) — K{ujbt + (wt + i^t) 



(28) 



IOCOS 



jet(^BT + lUT + U)\t) - j e t(^BT + Ujt) 

jet(uBT + Iujt) - j et (uJBT + Iujt - uj\t) 



(29) 



and 



K{ujbt + Iujt + uj\t) — K{u>bt + Iujt) 
K{ujbt + Iujt — uj\t) — K(ujbt + Iujt) 



(30) 



If the superlattice is unbiased (u>b — 0) and n is odd, we get for both frequency 
relations uj\ ± ui 2 = nui that -4™ h — A c °^ — 0. Thus the nondegenerate parametric 
amplification can be achieved in unbiased superlattice only if n is even. This re- 
sult is expected because the oscillations of the effective electron mass in unbiased 
superlattice take place only at even harmonics^!. 

3.4. The limits of weak and strong dissipation 

In the limit of weak dissipation r — > oo, we found that the incoherent absorption (|20p 
vanishes. On the other hand, by using the expressions for the coherent absorptions 
-4™ h and *4^° h ([2Tj) we can show that in this limit they satisfy the last equation in 
the set of Manley-Rowe relations (fT7| . In a similar way, we can also show that the 
coherent absorptions (|26|) lead to the the last equation in Manley-Rowe relations 

(LTHD- 



i 



i 
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In the limit of strong dissipation r — > (quasistatic limit), the situation is 
drastically different. Instead of behaving like a nonlinear reactance, the superlattice 
in this limit behaves like a nonlinear resistor. In the quasistatic limit, the asymptotic 
techniques (Appendix B in Ref. [22] and Ref. [14]) can be directly applied to find the 
formulas for the coherent and incoherent parts of absorptions from their general 
expressions @ The incoherent absorptions become 

^mcoh _ 1 ^ ^ (-Epump (t) I E cr ) } t , 

AZ° h = ^^-(3et(E„(t)/E cr )) t , (31) 

where E cr = h/edr is the Esaki-Tsu critical field corresponding to the peak current 
and j e t{x) is defined in We see that the incoherent absorptions are propor- 
tional to the dc differential conductivity, and thus in conditions where the domain 
formation is suppressed they are always positive!^ 

Similarly the coherent absorptions in the quasistatic limit can also be represented 
as integrals. In the case of u)\ + ll>2 = nui we get from formulas (f2T|) , (f22|) , (|23|) , (I24|) 
and (USD tnat 



Whereas in the case of w\ - lo 2 = ncu, formulas (J2BJ), ([27]). ([25]). ([21?]) and ||30J) are 
transformed to 



= ^^"^ ^{UE„{t)/E cr ) cos(n^)) t , 



We can make two important remarks concerning these equations. First, they are 
equivalent to the corresponding equations in the case of degenerate parametric 
amplification^^. Thus we can immediately say that the parametric gain can over- 
come the incoherent absorption leading to amplification of the probe fields in con- 
ditions where the formation of electric domains is suppressed. Secondly, we have 
an interesting relation for the coherent parts of the generated/absorbed power 
VZf = EiA c u f (i = 1,2). Namely 

Kt = Kf- (34) 

Finally, we would like to outline the differences in amplification/attenuation 
of the probe fields for the two mixing relations uj\ ± 0J2 — nu> in a variation of 
the dissipation strength. In the case w\ + 0J2 = nu>, the simultaneous amplification 



a Same results of course can be derived within the simple quasistatic theory! 
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of both modes is possible both in the limit of weak (|17p and strong dissipation 
([34]) . The main difference in this case is the change in the phase-dependence of 
absorptions from sine (fTS"]) . (fT6| to cosine (|32[) with an increase of dissipation. In 
the case u>i — L02 = nu), the changes are more drastic. Whereas the simultaneous 
amplification of both modes is possible in the limit of strong dissipation (I34p . it is 
forbidden in the limit of weak dissipation (fTS)) . 

4. Conclusion and perspectives 

We contributed to the theory of nondegenerate parametric amplification in super- 
lattices. We showed that the solution of Boltzmann transport equation for a tight- 
binding lattice satisfies the Manley-Rowe relation in the limit of weak dissipation. It 
demonstrates that we are really dealing with the parametric processes. Our theory 
predicts that in semiconductor superlattices the amplification of electric field can 
potentially be reached in a wide range of frequencies from microwaves up to high 
terahertz. 

We believe that similar parametric effects can be realized in various other physi- 
cal systems. Here we would like to mention only two of them: Nitride semiconductors 
and optical lattices. Since the dependence of effective mass o n the e nergy density 
in the lower part of the conduction band of dilute nitride 

alloysESEl 

resembles the 

corresponding dependence occurring in the miniband of superlattices, we speculate 
that our theories of nondegenerate and degenerate parametric amplification can be 
also applicable to these bulk semiconductors. On the other hand, the tight-binding 
optical lattices filled with cold atoms or Bose-Einstein condensates can demonstrate 
nonlinear tr anspor t in a single energy band in the both dissipative^^ and nondissi- 
pative Therefore the optical lattices can be considered as a playground 

for the demonstration of ac-driven transport satisfying the Manley-Rowe relations. 
A detailed study of the parametric resonance in these interesting systems goes be- 
yond the scope of the present paper. 
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